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Abstract  —  This  paper  presents  an  alternative  approach  to 
pulse  compression  that  is  based  upon  iterative  Minimum  Mean- 
Square  Error  (MMSE)  estimation.  It  is  similar  to  the  well- 
known  least  squares  (LS)  approach  but  does  not  suffer  from  the 
adverse  effects  caused  by  scatterers  closer  than  some  nominal 
range.  This  results  in  a  more  robust  estimate  of  the  radar 
returns  while  maintaining  nearly  the  same  sidelobe  level. 
Furthermore,  the  proposed  pulse  compression  scheme  is  highly 
parallelizable  in  range  and  therefore  can  be  computed 
efficiently. 


return  sample  can  be  expressed  as 

y(£)  =  xr(£)  s  +  v(£),  (2) 

where  x(f)  =  [x(f)  x(£  —  l)  x(f-./V  +  l)]7  consists  of 
samples  of  the  true  radar  channel  and  v(f)  is  additive  noise. 
The  matched  filter  output  can  therefore  be  written  as 

%W  =  sB[A(<)]Ts  +  s%(l),  (3) 


I.  Introduction 

Pulse  compression  allows  a  radar  to  obtain  the  range 
resolution  of  a  short  pulse  without  the  need  for  very  high 
peak  transmit  power.  This  is  accomplished  by  transmitting  a 
long  pulse  that  is  phase  or  frequency  modulated  to  generate 
a  wideband  signal.  The  wideband  signal  is  reflected  back  to 
the  radar  by  scatterers  in  the  beam  of  the  radar  which  can  be 
viewed  as  the  convolution  of  the  wideband  signal  with  a 
channel  that  is  representative  of  the  range  profile  illuminated 
by  the  radar.  The  purpose  of  pulse  compression  is  then  to 
estimate  the  radar  channel  based  on  the  known  transmitted 
signal  and  the  received  radar  return  signal. 

The  classical  approach  to  pulse  compression  is  known 
as  matched  filtering  [1]  which  has  been  shown  to  maximize 
the  received  signal-to-noise  ratio  (SNR)  and  is  accomplished 
by  convolving  the  transmitted  signal  with  the  received  radar 
return  signal.  One  can  represent  matched  filtering  in  the 
digital  domain  as  the  filtering  operation 

xMFU)  =  sHy(£),  (1) 

where  xMF(l) ,  for  (  =  0,-,L-l,  is  the  estimate  of  the  f1*1 
delayed  sample  (range  gate  index),  s  =  [.s,  s2  ■■■  s^]7  is  the 
length-A  sampled  version  of  the  wideband  transmitted 
waveform,  y(/j  =  [v(/j  y{£  + 1)  •••  y(/  + /V-l)]7  is  the 
length-N  vector  of  received  radar  return  samples,  and  (•) 
and  (»)7  are  the  complex  conjugate  (or  Hermitian)  and 
transpose  operations,  respectively.  Each  individual  radar 


where  v(f)=[v(0  v((  +  l)  •••  v(f  +  lV-l)]7  and 

x(()  x{(+\)  ■■■  x(t+N-\) 

x(7-l)  x(()  : 

:  x(l+ 1) 

_x{£-N+ 1)  •••  x(£-l)  x{t) 


[A(f)]  = 


(4) 


is  a  collection  of  sample-shifted  snapshots  of  the  radar 
channel. 

From  (4),  it  is  obvious  that  estimation  via  matched 
filtering  will  suffer  from  range  sidelobes  due  to  the  influence 
from  neighboring  radar  channel  coefficients  ( i.e .  the  off- 
diagonal  elements  of  [A(/')])-  To  alleviate  this  well-known 
effect,  least-squares  (LS)  solutions  have  been  proposed  [2]- 
[6]  that  decouple  neighboring  channel  coefficients  which 
have  been  smeared  by  the  transmitted  waveform.  The  LS 
solution  models  the  length-(L+A-l)  received  radar  return  as 

y  =  [S]x+v,  (5) 

where  x  =  [x(0)  x(l)  •••  x(L-l)]7  are  the  L  true  radar 
channel  coefficients,  v  =  [v(0)  v(l)  •  •  •  v(L  +  N  -  2)] 7  are 
additive  noise  samples,  and  the  convolution  of  the 
transmitted  waveform  with  the  radar  channel  is 
approximated  as  the  matrix  multiplication 
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The  general  form  of  the  LS  solution  is 


(6) 


xu=(s"S)-Vy. 


(7) 


For  the  received  signal  model  of  (5),  it  can  be  shown 
that  the  LS  solution  of  (7)  is  optimal  in  the  mean-square 
error  (MSE)  sense  when  the  additive  noise  is  white. 
However,  upon  further  inspection  one  finds  that  the  LS 
received  signal  model  does  not  completely  characterize  the 
received  radar  returns  because  it  does  not  account  for 
convolutional  components  involving  samples  of  x{£)  prior 
to  l  =  0  .  These  samples  fold  over  into  the  desired  window 
of  ranges  due  to  the  temporal  (and  hence  spatial)  extent  of 
the  transmitted  waveform.  The  result  is  that  the  presence  of 
a  significant  scatter  within  iV-1  range  gate  samples  prior  to 
4o)  can  have  a  significant  deleterious  effect  on  the 
estimation  of  a  large  portion  of  the  range  window. 

To  compensate  for  this  shortcoming  of  the  LS  estimator 
we  first  revisit  the  matched  filter  and  note  that  when  £  =  0  in 
(4),  the  matched  filter  does  account  for  the  samples  just  prior 
to  the  beginning  of  the  range  window.  It  is  therefore 
necessary  to  obtain  a  middle  ground  between  the  matched 
filter  and  the  LS  estimate.  The  remainder  of  this  paper 
develops  the  Minimum  Mean-Square  Error  (MMSE) 
estimator  [7]  for  pulse  compression  that  provides  a  robust 
estimate  of  the  radar  channel  while  maintaining  range 
sidelobes  at  or  very  near  the  level  attained  by  the  LS 
estimator. 


II.  MMSE  Pulse  Compression 

Minimum  Mean-Square  Error  estimation  is  a  Bayesian 
estimation  approach  that  employs  prior  information  in  order 
to  improve  estimation  accuracy.  The  exact  form  that  the 
prior  information  will  take  for  pulse  compression  will  be 
addressed  shortly.  First,  however,  the  signal  model  must  be 
constructed.  From  (3),  we  see  that  the  collection  of  N 
samples  of  the  received  radar  return  can  be  expressed  as 

y(f)  =  [A(f)]rs+v(f).  (8) 

This  is  the  same  received  signal  model  used  by  the  matched 
filter  and  takes  into  account  all  the  necessary  radar  channel 
components  for  a  given  received  return  sample.  To  develop 


the  MMSE  filter,  the  matched  filter  sH  in  (3)  is  replaced 
with  the  MMSE  filter,  denoted  v/H  (£) ,  in  which  the  exact 
form  of  the  MMSE  filter  is  dependent  upon  the  swath  of 
range  under  consideration.  Thereafter,  the  MMSE  cost 
function 


J(f)=£[|*OT-w"(f)yW|2]  (9) 

is  solved  for  each  range  gate  £  =  0,  •  •  • ,  L  —  1 ,  where  £[•] 
denotes  expectation.  This  is  done  by  differentiating  with 
respect  to  ww(f)  and  then  setting  the  result  to  zero.  The 
MMSE  filter  is  found  to  take  the  form 

w(f)=(£[y(f)  y'V)])'1^)  4(f)],  (10) 

where  (•)'  is  the  complex  conjugate  operation.  After 

substituting  for  y(()  from  (8)  and  assuming  that  the  radar 
channel  range  gates  are,  in  general,  uncorrelated  with  one 
another  and  are  also  uncorrelated  with  the  noise,  one  obtains 

w(f)  =  /7(f)([C(f)]  +  [R])'1s,  (11) 


where  p{£ )  =  E[  \x(£)f  ]  is  the  expected  power  of  x(l) ,  and 
[R]  =  E[\(£)  v"(f)]  is  the  noise  covariance  matrix.  On 
assuming  neighboring  range  cells  are  uncorrelated,  the 
(i,  j)th  element  of  the  matrix  [C(f)|  is 


[C  (£)]„j=E 


—  n  +  i  —  1)  s(n)  s"(n  —  i+  j) 


(12) 


in  which  kl  =  max{0,  i  —  j}  is  the  summation  lower  bound 
and  Kjj  =  minjiV  - 1,  N  - 1  +  i  -  j}  is  the  upper  bound.  Also, 
any  prior  information  regarding  the  noise  can  be  employed 
via  the  noise  covariance  matrix  [r]  .  For  instance,  for  a 
white  noise  assumption,  [R]  is  diagonal. 

Obviously,  in  its  current  state  the  MMSE  pulse 
compression  filter  is  a  function  of  the  powers  of  the  radar 
returns  in  surrounding  range  gates,  which  in  practice  are 
unavailable.  This  lack  of  prior  knowledge  can  be  taken  into 
account  by  setting  all  the  initial  radar  channel  estimates 
equal.  Therefore,  the  initial  MMSE  pulse  compression  filter 
reduces  to  the  form 


w  s[c]_1  s  (13) 

where  the  noise  term  is  assumed  negligible  and 


[c] 


Uj' 


X  s(n)  s*(n  -  i  +  j) 

n=rcL 


(14) 


is  invariant  to  the  range  gate  delay  i .  The  initial  MMSE 
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pulse  compression  filter  can  be  pre-computed  and  then 
implemented  in  the  same  way  as  the  traditional  matched 
filter.  The  inclusion  of  the  matrix  [c]  serves  to  provide  a 
“local”  LS  initial  guess. 

Figure  1  illustrates  the  MMSE  pulse  compression 
algorithm  for  three  iterations.  In  general,  the  reiterative 
MMSE  filter  operates  as  follows: 

1)  Collect  the  received  samples  y(-(M  — 1)( —  1)), — , 
y(L  +  M (N  -1)) ,  which  comprise  the  length-L  radar 
channel  window  along  with  the  (M  -  1  )(N  -  1)  samples 
prior  to  the  window  and  the  M(N  -  1)  samples  after  the 
window. 

2)  Apply  the  initial  MMSE  pulse  compression  filter 
from  (13)  to  obtain  the  initial  radar  channel  estimates 
^  {-(M  -  2)(N  - 1) ),  ■  •  • ,  x,  (L +  (M  -  1){N  - 1)) . 

3)  Compute  the  initial  power  estimates  pi(/)  =  \xl(i)'\ 
for  l  =  -(M-2)(  N- 1),-  -,L  +  (M-1)(JV-1)  which 
are  used  to  compute  the  range-dependent  filters  W[(f) 
as  in  ( 1 1),  and  then  apply  to  y(l)  to  obtain  x2{i)  ■ 

4)  Repeat  2)  and  3),  changing  the  indices  where 
appropriate,  until  the  desired  length-L  range  window  is 
reached. 


Fig.  1.  Operation  of  reiterative  MMSE  pulse  compression  filter 


The  initial  guess  of  the  radar  channel  found  by  applying 
the  MMSE  filter  can  be  used  as  a  priori  information  to 
reiterate  the  MMSE  filter  and  improve  performance.  This  is 
done  by  employing  the  MMSE  filter  formulation  from  (11) 
in  which  the  respective  powers  of  the  radar  returns  are  taken 
from  the  current  estimate.  It  has  been  found  that  two  or 
three  reiteration  steps  allow  the  MMSE  pulse  compression 
filter  to  exclude  the  effects  of  scatterers  prior  to  the  range 
window,  as  well  as  to  suppress  the  range  sidelobes  very 


close  to  the  level  of  the  noise  floor.  The  MMSE  filter  does 
especially  well  when  the  radar  channel  is  somewhat  sparsely 
parameterized  (i.e.  highly  spiky),  as  is  the  case  with  high 
range  resolution  radar. 

Note  that  each  reiteration  step  will  reduce  the  number  of 
range  gate  estimates  by  2(lV-l).  To  counteract  this,  it  is 
necessary  to  increase  the  size  of  the  range  window  by 
2M(N  -l)  samples,  where  M  is  the  number  of  reiteration 
steps.  Typically,  however,  L»  N  so  that  this  reduction  in 
range  window  size  is  negligible. 


III.  Implementation 

One  may  immediately  notice  that  in  order  to  reiterate 
the  MMSE  pulse  compression  filter  it  is  necessary  to  invert 
an  N  x  N  matrix  for  every  range  gate  which  requires  on  the 
order  of  N 3  operations  for  each  inverse.  This  can  cause  a 
computational  burden  whenever  either  N  or  L  is  large. 
However,  the  inversion  of  the  roughly  L  matrices  can  be 
divided  into  small  groups  and  computed  in  parallel.  For  this 
reason,  the  computational  complexity  of  the  MMSE  filter 
scales  with  the  number  of  parallel  processors  and  hence  can 
be  computationally  feasible  given  enough  processors. 

Another  important  factor  regarding  practical 
implementation  is  the  non-singularity  of  the  N  x  IV  matrix 
([CCO]+[R]).  This  can  be  addressed  by  instituting  a 
nominal  level  for  which  the  estimated  returns  are  not 
allowed  to  go  below.  An  alternative  to  this  would  be  to  re- 
estimate  only  those  range  gates  that  are  above  some 
threshold  since  the  small-valued  range  gates  most  likely  do 
not  contain  a  target  and  are  therefore  of  little  interest. 


IV.  Simulation  Results 

To  demonstrate  the  performance  of  the  MMSE  pulse 
compression  filter  we  examine  two  cases.  The  first  case  is 
typical  of  the  scenario  often  addressed  for  LS  estimation 
techniques  and  consists  of  a  radar  channel  with  a  single  large 
target  return  in  noise.  The  other  case  we  address  is  when 
there  is  a  second  large  target  that  resides  in  a  range  just  prior 
to  the  range  window.  In  this  case  it  is  expected  that  the  LS 
estimator  will  substantially  degrade  since  this  important 
region  is  not  accounted  for  in  the  received  signal  model. 
The  waveform  selected  is  the  length  N  =  30  polyphase 
modulated  Lewis-Kretschmer  P3  code  [8],  which  on  receive 
(after  down-conversion  to  baseband)  is  defined  as 


s{ri)  =  exp 


n  =  0,1,  •  •  • ,  V  - 1 .  (15) 


The  noise  and  clutter  are  modeled  as  zero-mean  Gaussian 
with  powers  set  to  -50  dB  and  -30  dB,  respectively.  For 
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both  cases  we  perform  two  iterations  of  the  MMSE  filter  and 
then  compare  the  MMSE-estimated  radar  channel  with  the 
true  radar  channel,  as  well  as  with  the  results  obtained  from 
using  LS  and  the  matched  filter. 

The  radar  channel  for  the  first  case  is  depicted  in  Fig.  2 
for  a  single  target  present  in  the  range  window  which 
consists  of  200  range  gates.  Figure  3  illustrates  the  results 
from  the  different  estimation  techniques  in  which  the 
matched  filter  experiences  significant  range  sidelobes  while 
both  LS  and  MMSE  have  suppressed  the  range  sidelobes  so 
that  the  true  target  location  is  evident. 


range  gate  index 


Fig.  2.  Radar  channel  with  single  target  in  the  range  window 


range  gate  index 


Fig.  3.  Performance  of  MMSE,  LS,  &  matched  filter 
for  single-target  channel 

For  the  scenario  just  described,  the  LS  and  MMSE 
estimators  perform  almost  identically.  Flowever,  when  there 


performance  of  the  MMSE  estimator  is  indistinguishable 
from  the  previous  case  where  only  a  single  target  was 
present.  This  is  obviously  due  to  the  fact  that  the  MMSE 
estimator  takes  all  necessary  range  regions  into  account 
when  estimating  the  radar  channel. 


range  gate  index 


Fig.  4.  Radar  channel  with  one  target  in  range  window  and 
one  target  just  before  range  window 


range  gate  index 


Fig.  5.  Performance  of  MMSE,  LS,  &  matched  filter 


for  target  before  range  window 


Table  1  presents  the  Mean  Squared-Error  (MSE)  for 
both  radar  channels  discussed  using  LS,  MMSE,  and 
normalized  matched  filtering.  The  MSE  is  averaged  for  all 
200  range  gates  and  over  100  runs  using  the  same  target(s), 
with  the  clutter  and  noise  distributed  according  to  a  zero- 
mean  Gaussian  distribution  for  each  run  independently. 


is  a  significant  scatterer  present  just  prior  to  the  range 
window,  as  depicted  in  Fig.  4,  the  LS  estimator  is  expected 
to  degrade  substantially.  From  Fig.  5,  the  LS  estimator  truly 
does  suffer  from  severe  mis-estimation  when  a  significant 
scatterer  cannot  be  expressed  in  the  model.  However,  the 


Table  1.  MSE  performance  comparison 


— _____ 

Norm.  MF 

LS 

MMSE 

Channel  1 

-12.3  dB 

-40.9  dB 

-35.8  dB 

Channel  2 

-11.2  dB 

-11.3  dB 

-34.7  dB 
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For  channel  1  (single  target),  the  normalized  matched 
filter  has  a  large  MSE  due  to  range  sidelobes,  while  LS  and 
MMSE  perform  nearly  the  same  with  LS  just  marginally 
better.  However,  for  channel  2  (target  just  prior  to  the  range 
window),  the  MSE  attained  by  LS  degrades  to  nearly  that  of 
the  normalized  matched  filter,  while  the  MMSE  maintains 
roughly  the  same  MSE  as  in  the  previous  case. 

V.  Conclusion 

Least-Squares  pulse-compression  estimation  of  the  radar 
channel  is  a  well-known  approach  that  has  been  shown  to  be 
an  optimal  estimator  for  the  assumed  radar  channel  model. 
However,  scatterers  may  exist  prior  to  the  range  window  of 
interest  which  would  inherently  invalidate  the  LS  radar 
channel  model.  To  alleviate  this  problem,  the  MMSE 
estimator  is  proposed  that  reiteratively  finds  a  “local"  least- 
squares  solution  and  in  so  doing  takes  all  necessary  radar 
channel  information  into  account.  The  result  is  an  estimator 
that  substantially  reduces  range  sidelobes  compared  to 
matched  filtering,  as  well  as  provides  robustness  to  scatterers 
that  may  lie  outside  the  range  window  of  interest. 
Furthermore,  reiterative  MMSE  estimation  is  not  limited  to 
pulse  compression  alone.  It  may  also  find  use  in  range 
profiling,  image  recognition  for  SAR  and  ISAR,  or  any  other 
application  requiring  robust  deconvolution  of  a  known 
waveform  with  an  unknown  impulse  response. 
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